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I problem. Our work integrates and extends previous work by Gudkov et al. (,cond-mat/0209112 ) and 
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5h . We investigate classical and quantum physics-based algorithms for solving the graph isomorphism 
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by Rudolph (quant-ph/0206068 |. Gudkov et al. propose an algorithm intended to solve the graph 



■ isomorphism problem in polynomial time by mimicking a classical dynamical many-particle process. 

> : 

, We show that this algorithm fails to distinguish pairs of non-isomorphic strongly regular graphs, 



thus providing an infinite class of counterexamples. We also show that the simplest quantum 
generalization of the algorithm also fails. However, by combining Gudkov et al.'s algorithm with 
a construction proposed by Rudoph in which one examines a graph describing the dynamics of 



o 

Oh, two particles on the original graph, we find an algorithm that successfully distinguishes all pairs 



' of non-isomorphic strongly regular graphs that we tested (with up to 29 vertices). 
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Introduction. The graph isomorphism problem plays a central role in the theory of 
computational complexity and has importance in physics and chemistry as well^, 0]. A 
graph is a set of points, or vertices, and a set of edges, or unordered pairs of those 
points. If two graphs differ only in the labelling of their points, then we say they are 
isomorphic, otherwise they are non-isomorphic. The graph isomorphism problem is to 
determine whether or not there is an algorithm that runs in polynomial time {t ~ A^^, with 
X independent of A^) that distinguishes non-isomorphic pairs with certainty. The efficient 
enumeration of possible distinct atomic clusters of size A^, in which a cluster is defined by 
the bonds between its atoms, is the same problem in another guise 0]. 

Typical instances of graph isomorphism (Gl) can be solved in polynomial time be- 
cause two randomly chosen graphs with identical numbers of vertices and edges typically 
have different degree and eigenvalue distributions. Moreover, GI can be solved efficiently 
for restricted classes of graphs, such as trees planar graphs yi, graphs with bounded 

nn n q 

degree p, 131, bounded eigenvalue multiplicity [8|], and bounded average genus |9|]. However, no 
algorithm for solving GI for all graphs is presently known; the best existing upper bound is 
exp There is evidence that indicates that GI is not NP-complete: first, count- 



10, 



ing the number of isomorphisms is reducible to the decisional version of the problem un- 
like the situation for all known NP-complete problems, and second, if GI were NP-complete 
then certain complexity classes that are believed to be different would be identical^, Q- 
Thus, it is believed that either GI is in P or else it is in the class of problems that are neither 

n 

in P nor are NP-complete [1^ 

One way to solve GI is to solve the hidden subgroup problem for the permutation 
group. Unfortunately, though the hidden subgroup problem for abelian groups can be 
solved efficiently^^, no efficient alg orithm for solving the hidden subgroup problem for 



the permutation group is known 



16 



In this paper we investigate classical and quantum approaches to solving the graph iso- 
morphism problem that are motivated by physical processes, jisl Our work combines and 
extends ideas in Gudkov et a/. [22] and by Rudolph j22| for attacking the GI problem using 
algorithms based on physical intuition. We show that the interesting classical dynamical 
algorithm proposed by Gudkov et al. fails to distinguish an infinite set of pairs of non- 
isomorphic graphs, and thus does not solve the GI problem in polynomial time. We trace 
this failure to certain algebraic properties of these particular pairs of graphs and show that 



the simplest quantum generalizations of the Gudkov et al. algorithm also fails to distinguish 
these pairs of graphs. However, an algorithm obtained by combining the basic idea of the 



Gudkov et al. algorithm with a construction proposed by Rudolph j23| does distinguish all 
pairs of graphs that we have examined, including many with the same eigenvalue spectra. 

The detailed statement of the GI problem is as follows. We are given two graphs. The 
first is a set of vertices {fi,f2, ■■■Vn} , together with a set of edges, or unordered pairs 
{ei,e2,...} connecting pairs of these vertices, while the second graph is a set of vertices 
{v[,V2, --.v'^} , and a set of edges {e[,e'2, ...}. Each is associated with a pair {va,Vb}, 
and each e[ is associated with a pair {v'^,v'f^} . We wish to determine whether there exists 
a permutation V of the Vi such that the set of pairs {Vva,Vvb} is identical with the set of 
pairs {v'^,v'i,}. 

An equivalent but more useful formulation for our purposes is to represent each graph 
by its adjacency matrix A. Aab is an N ^ N matrix such that Aab = 1 if there is a pair 
e = {va,Vb} and Aab = otherwise. Two graphs are isomorphic if and only if there exists 
a permutation matrix P such that A' = PAP^^. 

Classical algorithm. In the Gudkov et al. approach to this problemj2^, each 
vertex of the graph is mapped into a point in an A^— dimensional vector space. The vertex 
fi sits initially at the point fi = (1, 0, 0, 0), V2 at r2 = (0, 1, 0, ...0), and so on. These are 
the vertices of an (A^ — 1)— dimensional simplex. We now view these as mass-points with 
forces acting pairwise among them. There is one force law if two particles are connected by 
an edge and a different force law if they are not. The forces are derived from potentials Ui 
and U2 that depend only on the distances between the points in the pair. From the initial 
configuration, the particles move in time according to the relaxational equations of motion 

where 

^a = -V^„f/(ri,F2,...) 

and 

U = J2 AabU,{\fa - n\) + - Aab)U2{\fa - (l) 

a>b a>b 

The motion may be computed by any convenient algorithm. After a time T the positions 
are given by ra{T) for a = 1,2,...N. We now compute the set of A^(A^ — l)/2 distances 



dab — \^a{T) — rb{T) \ . Given a second graph, we compute = l^aC^) ~ using the 

same prescription. The sets {dab} and {rf^b}) being non-negative real numbers, may be 
arranged in increasing order. If the resulting sets are identical, then it is conjectured that 
the graphs are isomorphic. Note that as long as T is not an exponentially long time, 
this is a polynomial-time algorithm. The sets dab and rf^^ can be computed, ordered, and 
compared efficiently. We wish to examine the claim that if {dab} and {d'^f,} are the same 
up to reordering, then their parent graphs are isomorphic. 

For purposes of clarity, we shall initially consider a slightly simpler model, one in which 
a harmonic force acts only between particles connected by edges: 

U = -fiJ^'^^bK-fblV^. 

a>b 

We then define an TV x TV matrix X, where the i-th coordinate of the a-th "particle" is 
denoted by Xai- Thus the above initial condition can be rewritten as 

Xai{t = 0) = 6ai, 

where 5ai is the Kronecker symbol: 5ai — 1 H a — i and 5ai — otherwise. The equations 
of motion are 

, = ^ -^ife'' = ^ Aab{Xai — Xbi) — — AabXbi + ^ ^ Aac5abXbi 
b b b b c 

— LgbXbi, 

b 

where 

Lab = Sab ^ ^ ^ac ~ ^a6 = Dab ~ ^ab, 

c 

is the Laplacian matrix. The diagonal matrix D is the degree sequence matrix: Daa is the 
number of edges incident to the vertex a. Note that L is symmetric: L = . The algebraic 
isomorphism criterion mentioned in the introduction may also be put in terms of L : two 
graphs defined by the Laplacian matrices L and L' are isomorphic if and only if there exists 
a permutation matrix P such that 

L' = PLP^. 

The solution to the equation of motion for the particles defined by U is 

Xa^{t) = ie^')^^Xb.{t = 0) = {e^^)^^ 



so the final positions are 

X(T) = e^'^, 

in a matrix notation. In order to compute the distances, we note that the set of dot products 
between the position vectors may be written as: 

Sab = ^ai^bi 
i 

which is the matrix 

S = XX^""^P°"" = X^ = e^^* = l + 2tL + ^^L^ + ... 
Another graph would be characterized by a different dot product matrix 

S' = e^^'*. 

Since the N"^ distances dab satisfy d^f^ = Saa + Sbb — 2Sab, comparing the set of numbers in the 
matrices S and S' is essentially the same as comparing the distances. (We shall consider 
the relationship more carefully below.) 

If the graphs are isomorphic, then clearly S' is a rearrangement of S : 

S' = e^^'^ = e^PL""-'^ = Pe'L'p-' = PSP-\ 

The interesting question is whether the converse is also true. 

In fact there do exist many interesting graph pairs for which the algorithm works. Because 
S is most easily computed by diagonalizing the real symmetric matrix L, it is natural to ask 
whether pairs of non-isomorphic but isospectral graphs can be distinguished by the method. 
Isospectral graphs are those for which the eigenvalues of A and A' are the same. We have 
investigated this question for some small graphs that are isospectral but not isomorphic. A 
simple illustrative pair is shown in Fig. 

We take the total time interval as T = 1 and compute numerically the dynamics of the 
simple harmonic model for 10 steps of length 0.1 using the first-order Euler algorithm and 
finally obtain a normalized X{T). The sorted (i^^ and d^'j, for two graphs G and G' in Fig. ^ 
are respectively 

dlb = {[0]^[0.0785]l^[3.9685]«}. 

ci2; = {[-1.9216]^[-0.2570]^[0]^[0.1406]«}. 



G G' 
FIG. 1: Two isospectral graphs 
In these expressions the superscripts denote the multiphcity of the number in square brackets. 



Thus the algorithm of Ref. [2^ works for this non-isomorphic isospectral pair. The 
dynamical algorithm also distinguishes successfully some pairs of graphs that have both 
identical degree distributions and identical Laplacian spectra[29]. 

It has been known for decades, however, that certain classes of graphs are difficult to dis- 
tinguish by elementary methods. An important intransigent class is the so-called "strongly 
regular graphs" (SRG's)j24, A SRG with parameters {N,k,\,fi) is a graph with 
vertices in which each vertex has k neighbors, each pair of adjacent vertices has A neighbors 
in common, and each pair of non- adjacent vertices has /i neighbors in common. An example 
known as L2(3) with parameters (9, 4, 1, 2) is shown in Fig.|21 Many pairs of nonisomorphic 
SRGs with the same parameter set are known. 30] 

The adjacency matrix A of a SRG has some interesting algebraic properties. For a 
general graph, the (a, b) entry of is the number of vertices adjacent to both a and b. For 



SRGs, this number is {A^ 



lab 



kii a = b, {A^ 



'ab 



A if a is adjacent to b, and {A^ 



'ab 



/i if 



a is not adjacent to b. Hence 



A^ = kl + \A + fi{J - I - A), 
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FIG. 2: L2(3) (9,4,1,2) 



where / is the identity matrix and J is the matrix consisting entirely of I's. = NJ. A 
and J also have the properties that 

AJ^JA^ kJ. 

This follows from the fact that multiphcation of a matrix by J has the effect of adding the 
rows or columns of the matrix. For A, this sum is just the number of neighbors. The 
Laplacian matrix for a SRG is 

L^kl -A 



These equations show that the three matrices {/, J, L} form a commutative, associative 
algebra. The elements of the algebra have the form R — fl + gj + hL, where /, g and h 
are real numbers. The multiplication rule is: 

RE! = E!R = (// + ^ J + hL) {f'l + g'J + h'L) 

= [ff - [k^ -k{\-^i + l)+ ii]hl^] I + ifg' + gf + Ngg' + i,hh')J 
+ [fh' + hf + {2k + fi- \)hh'\ L 



The structure of the algebra is therefore independent of the precise form of the L matrix, 
depending only on the {N, k, A, /i) parameters. 

We now consider two non-isomorphic graphs characterized by Laplacian matrices L and L' 
that share a set {N, k, A, //) . Let the corresponding dot product matrices be 5" = exp(2LT) 
and S' = exp(2L'T). Since L and L' define the same algebra, we have 

S {N, k, A, i^,T) = f {N, k, \,iJ,T)I + g {N, k, X,n,T)J + h {N, k, A, /x, T) L 

and 

S' (N, k, A, ii,T)^ f (N, k, X,ii,T)I + g (N, k, X,ii,T)J + h (N, k, A, T) L' 



where the /, g, and h are definite functions of the parameters — the point being that the 



functions are the same for the two graphs. More exphcitly, we have that 

/ f + g + kh g-hAu g - hA^ \ 
g - hA2i f + g + kh g - /lAas 
g - hAsi g - hA32 f + g + kh . . . . 

S 



\ ■ / 

The same relation holds when S is replaced by S' and A is replaced by A'. In any row or 
column of A or A', precisely k entries are equal to 1 and N — k entries are equal to 0. It 
now follows that, considered as a set of numbers, 5* has entries equal to f + g + kh, Nk 
entries equal to g — h, and A^(A^ — A; — 1) entries equal to g. The same holds true for S'. 



The N"^ squared distances satisfy (d^ 



'ab 



Saa + Sbb - '2,8 ab = 2{f + kk) + 2hAab, and 



ab 



0. 



id%b = S'aa + S'bb - 2S'^, = 2{f + kh) + 2hA'^, with a ^ 6. For a = 6, {<P)^^ 
There are Nk nonzero entries of Aab and A'^^ with a ^ h. Hence, for both graphs there will 
be Nk distances 2(/ + kh) + 2h, N{N — k — 1) distances equal to 2{f + kh), and distances 
equal to 0. Thus it is impossible to distinguish this pair of non-isomorphic graphs by this 
simplified algorithm. 

To extend this to the actual algorithm of Gudkov et al, we examine their general equa- 
tions of motion ((T)). From a SRG defined by L, these equations produce a matrix X(T) 
which starts life as X(t = 0) = /. Any numerical solution of the differential equation is 
simply a finite sequence of matrix multiplications and additions. All such operations belong 
to the algebra defined by /, J, L. The dot product matrix S also belongs to this algebra. 
Thus the result is again characterized by only three numbers /, g, h. For a different SRG 
defined by L' but with the same parameter set, the algorithm produces the same f,g,h, 
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FIG. 3: Latin squares L-i{A) (16,9,4,6) 



26] 



since the algebra for the two graphs is the same. The above argument for the simphfied 
algorithm then goes through without further modification. 

To illustrate the breakdown of the algorithm, we shall apply it to the smallest pair of 
non- isomorphic SRGs. These are the "Latin square" graphs with size = 16. Latin squares 
are two-dimensional M x M arrays of the numbers 1 to M, arranged so that in each row and 
column no number is repeated. Two examples are shown in Fig. IHl Latin square graphs are 
constructed from Latin squares as follows: Given a Latin square of order M, the vertices 
are the N = cells. Two vertices are adjacent if they lie in the same row or column or 
if they share the same integer label. 

We use the following non-harmonic potential Ui to calculate normalized X{T) and (i^^ 
again using the first-order Euler algorithm for the two non-isomorphic Latin square graphs 
drawn from Fig. IHl 

f/i = -A J2ira - n)' + B Y,{ra - nf , 

where the sum is over pairs of connected vertices. Again taking time interval T = 1 with 
step of length 0.1, we obtain the distances d^f^ for A = 1, B = 1 

dlb = {[0]'^[1.8641]9^ [2.3129]^^^}. 
= {[0]^^[1.8641]9^ [2.3129]^^^}. 

f^'\ h^'^ are computed from c?^^'^ 

f = f = -1.0876 
h = h' = 0.2244 

g = g', but these quantities do not affect the distances. We have also verified that the two 
sets are identical at each of the discrete time steps. As one would expect, the multiplicity 
of each distinct distance depends only on A^ and k and are independent of time. 

For completeness we tried a different non-harmonic attractive(repulsive) potential Ui{fa— 
H)iU2{fa — r^)) whose force is expressed as 



Using this potential and still taking T = 1, we obtain the distances 

dl, = {[0]l^[1.4991]^^ [2.4486] 1^^}. 
dab = {[0]'^[1.4991]9^ [2.4486]^^^}. 

with /W, /iW 

f = f = -3.5482 
h = h' = 0.4748 

Thus, using a non-harmonic potential does not enable the dynamical algorithm to distinguish 
these graphs. 

Because we have identified two non-isomorphic grap hs that result in the same list of 



221. There are an infinite number of 



distances, we have disproved the conjecture of Ref. 
nonisomorphic pairs of SRGs with identical parameter sets, so the number of counterexample 
pairs is infinite. 

Quantum algorithms. The same argument can also be used to show that a simple 
one-particle quantum random walk algorithm also fails to solve the GI problem. Consider 
the vertices of the graph as states |j) in a Hilbert space. The Hamiltonian for the walk is 

H = -^Aabclcb, 

ab 

where the operator clcb is defined by (i |cjjCfe| j) = Sia^bj- In physics terms, this is simply a 
tight-binding model with bonds on the vertices of the graph. We now consider possible 
starting wavefunctions 

|^i(t = 0)) = |l), = 0)) = |2) , ... 

and we evolve these forward in time according to the usual time-dependent Schroedinger 
equation 

'^--W (2) 
for a time T. We then compute the N x N matrix of overlaps 

O., = (^,(0)|^,.(T)>. 

We might then conjecture that this matrix, considered as a set of complex numbers, is 
different for non-isomorphic graphs. It can distinguish isospectral graphs, since it uses 



information about the eigenvectors of iJ, not just information about the spectrum. However, 
our algebraic argument is trivially extended to the now complex algebra defined by /, J , L. 
For non-isomorphic SRG's with identical parameter sets, we again find that the two matrices 
of overlaps are the same after a rearrangement, so the conjecture is invalid. 



Though the two strong, 
by the algorithm of Ref. 



regular Latin square graphs of Fig. |H1 are not distinguished 
I, they can be distinguished in polynomial time by using a 



construction proposed by Rudolph j23|, in which the original adjacency matrix is interpreted 
as the Hamiltonian of a tight-binding model. The original graph with vertices describes 
the possible states (positions) of a single particle. However, we can also consider the 
quantum-mechanical motion of many particles on the same graph. Rudolph uses symmetric 
wavefunctions but forbids double occupancy of any site, corresponding to a hard-core boson 
model. Rudolph showed that the spectra of the 3-particle matrices obtained from two non- 
isomorphic regular graphs with identical single-particle spectra are different, demonstrating 
that the multiparticle construction does increase the power of the algorithm to distinguish 
similar graphs. 

Here we combine Rudolph's multiparticle construction with the dynamical algorithm 
for wavefunction overlaps; this hybrid algorithm has the advantage that it can distinguish 
nonisomorphic strongly regular graphs using the 2-particle matrices, as opposed to 3-particle 
matrices needed if the matrix eigenvalues are examined. In addition to Rudolph's original 
case of hard-core bosons, we also examine non-hard-core bosons and noninteracting spinless 
fermions. 

I — I 

For bosons we consider a simple Hubbard Hamiltonian [27], in which each boson can hop 
between two vertices if and only if the vertices are connected by an edge, and in addition 
there is an energy cost U if two bosons are on the same site. Using as a basis the states 
with particles at vertices i and j, the matrix K for the level 2 (i.e. two-particle) graph for 
bosons is specified by the Hamiltonian matrix elements 

Kl ^ - (zj \H\ kl) 



6uAkj + 6jkAii + 6ikAji + 6jiAik if ^ 7^ j and k ^ I 

U 6ik if i = j and k = I 

-^{6iiAkj + SjkAii + SikAji + SjiAik) if i = j or k = I, 
v2 



where for bosonic statistics we have A^(A^ + l)/2 initial two-particle states 

= 111) , |12) , |13) , ... \1N) , |22) , |23) , |24) \NN) . 

In the hard core limit U oc the basis states with doubly-occupied sites can be ignored, 
the Hilbert space has N{N — l)/2 dimensions, and the Hamiltonian matrix elements are 

Kg^,"! ^ - (.J \H\ kl) 

= SiiAkj + SjkAii + SikAji + SjiAik , 

where now we require i ^ j and k ^ I. 

For Fermi statistics, since two fermions cannot occupy the same vertex, the Hilbert space 
has iV(A^ — 1)/2 dimensions, and we can choose the basis states \ij) , i ^ j- The Hamiltonian 
matrix elements are 

^fj, kl = ^il^kj + ^jk^il — ^ikAjl — SjiAik- 

Thus we extend an adjacency matrix of rank to a matrix K of higher rank, either N{N+ 
l)/2 (for non-hard-core bosons) or A^(A^ — l)/2 (for hard-core bosons and for fermions). 
Technically, except for hard core bosons, the matrix K is not an adjacency matrix since it 
has elements other than and 1. For fermionic statistics, some entries of K are equal to 
— 1, while for soft-core bosons some entries are equal to -\/2 and others to U. K can be 
pictured basically as a matrix in which every off-diagonal element represents the probability 
amplitude for particles to hop from one state to another state. Accordingly, the information 
of the adjacency matrix A is embedded in the corresponding K. 

Again we evolve forward the initial two-particle states in time according to the quantum 
mechanical evolution Eq. Q and ask whether any pair of nonisomorphic graphs has two 
distinct sets of Oij{T). The test of isomorphism is whether the sets Oij and O-^- are the same 
after rearranging. It is of course possible to order these sets first by their real parts and 
then by their imaginary parts. A simple way to compare the ordered sets is to compute the 
numbers R and I defined as 

i?(T) = 5^ |ReOy(T)-Re6:j(T)| (3) 



/(T) = 5^|Im6ij(T)-Im6:j(T)| (4) 

id 



R, I 
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FIG. 4: Variation of the numbers R and I (defined in Eqs. |31and0]) as a function U (potential) for 
the two non-isomorphic Latin square graphs with N=16. As U goes to zero, R and I vanish. 

where Oij, 0[j are the elements of Oj^, O-^- ordered by their real parts, while Ojj, 0[j are 
ordered by their imaginary parts. If either R{T) or /(T) is nonzero, the graphs are non- 
isomorphic. The numbers R and / are sufficient to distinguish the non-isomorphic graphs 
in this paper. 

To probe this approach, we again take pairs of non-isomorphic strongly regular graphs, 
two examples being the pairs of non-isomorphic Latin square graphs in Figs. El and For 
each graph in a given pair, we find numerically the eigenvectors and eigenvalues of K and 
use them to calculate the Oij, O-^ and then R{T) and /(T) for T = 1. The qualitative 
behavior does not depend on the choice of T. Table IJ shows that for all these pairs of graphs 
both R and / vanish for noninteracting bosons, but both R and / are nonzero for hard core 
bosons and for noninteracting fermions. We find also that R and / are nonzero for graphs 
with nonzero but finite U . An illustration of this is shown in Fig. |3]the pair with N = 16. 
Note that the linear term in U is sufficient to distinguish these graphs. 

Thus, we find that all the non-isomorphic pairs of strongly regular graphs that we have 
examined (with size N = 16, 25, 26, 28, 29) can be distinguished by using two-fermion 
and the interacting boson algorithms. It is necessary to examine all matrix elements of K, 
not just the spectrum of K. The K matrices of the two non-isomorphic graphs still share 
the same set of eigenvalues. In physical terms, this is due to the additivity of energies for 




graph specification 


fermions 


noninteracting bosons 


hard-core bosons 


(16,9,4,6) 


R=1.37 


R=0 


R=110.66 




1=3.01 


1=0 


1=81.53 


(25,12,5,6) 


R=1.24 


R=0 


R=129.66 




1=1.93 


1=0 


1=198.53 


(26,10,3,4) 


R=1.91 


R=0 


R=14.88 




1=0.65 


1=0 


1=23.47 


(28,12,6,4) 


R=1.82 


R=0 


R=87.27 




1=1.25 


1=0 


1=95.11 


(29,14,6,7) 


R=3.50 


R=0 


R=28.69 




1=3.73 


1=0 


1=42.51 



TABLE I: Table of results for hybrid dynamical algorithm for pairs of nonisomorphic strongly 
regular graphs for non-interacting fermions, non-interacting bosons, and hard-core bosons. Each 
pair of graphs has the same parameters (A^, k, A, /x), where N is the number of vertices, each vertex 
has k neighbors, each pair of adjacent vertices has A neighbors in common, and each pair of non- 
adjacent vertices has fi neighbors in common. The algorithm distinguishes the nonisomorphic pairs 
when using either hard-core bosons or fermions, but not when using noninteracting bosons. 

non-interacting particles. 

Conclusion. Our work suggests two conjectures for the solution of the GI problem. 
The first, which we doubt but have not been able to disprove, is that the two-particle 
quantum algorithm with computation of overlaps solves GI. A more plausible conjecture 
is that the motion of Np quantum-mechanical fermions or interacting bosons on the graph, 
where Np ~ A^, say Np = N/2, solves GI. From the standpoint of efficiency, we must 
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FIG. 5: Latin squares L3(5) (25,12,5,6) 



distinguish whether the algorithm is run on a classical or a quantum computer. The two- 
particle algorithms already run in polynomial time on a classical computer. The Oij matrix 
contains 0{N'^) entries, each of which is computed in a time of 0(A^^). If this algorithm is 
sufficient to distinguish all graphs, then GI is in P. If, on the other hand, we must consider 
the motion of Np particles, where Np ~ N, then the time is roughly {^) ^ {n/2) ' "^hich is 
exponential in A^. The question of efficiency on a quantum computer is more interesting. 
The motion of Np N particles can be mapped onto the Heisenberg model of a spin 
system at a fixed magnetization, which is nothing more than the dynamics of qubits with 
a constraint. Thus the evolution takes place in polynomial time on a quantum computer. 
On the other hand, we need an exponential number of overlaps if a naive algorithm is chosen. 
However, since all we need to do is ask whether the evolution of the two systems is in some 
sense similar, there may exist a preparation that entangles the graphs and a measurement 
that captures the needed information. 

In conclusion, we have shown that strongly regular graphs provide a useful testbed for 
both classical and quantum dynamical algorithms that aim to solve the graph isomorphism 
problem. Pairs of graphs in this class cannot be distinguished by the classical algorithm of 
Gudkov et al. j^j, but quantum algorithms combi ning the dynamical evolution of Gudkov 
et al.'s algorithm with a construction of Rudolph [2j| using interacting bosons as well as 
noninteracting fermions can distinguish pairs of SRGs of order < 29, while a two-boson 

n 

noninteracting algorithm fails. 28] Interesting open questions include whether either a Np- 
fermion or A'p-hard-core-boson algorithm solves the GI problem, and if so, what value of 
Np suffices, and whether an algorithm with Np N can be implemented efficiently on a 
quantum computer. 

We thank Eric Bach and Dieter van Melkebeek for extremely useful conversations. This 
research was supported by the NSF QuBIC program under Grant No. NSF-ITR-0130400. 
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